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Measurements of sin2/3 and cos2/? from b — > ccs decays at BABAR. 
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Recent measurements of sin2/3 and cos2/3 using b — > ccs decays are presented using data collected 
by the BABAR experiment at the PEP-II asymmetric-energy B-factory. 



I. INTRODUCTION. 



II. sin2/3 FROM B° -> CHARMONIUM + K° 
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The Standard Model (SM) of particle physics de- 
scribes charge conjugation-parity (CP) violation as a 
consequence of a complex phase in the three-generation 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix [1]. In this framework, measurements of CP asym- 
metries in the proper-time distribution of neutral B de- 
cays to CP eigenstates containing a charmonium and K° 
meson provide a direct measurement of sin2/3 [2], The 
unitarity triangle angle (3 is arg [ — V r cd V^.* 3 /V td V t * J ] where 
the Vij are CKM matrix elements. CP violation in the 
B-meson system was established by the BABAR [3] and 
Belle [4] collaborations in 2001. 

The BABAR detector [5] is located at the SLAC PEP- 
II e + e~ asymmetric energy B-factory [6] where data is 
collected on or just below the T(4S) resonance. A small 
fraction (ps 10%) is collected at approximately 40 MeV 
below the T(AS) resonance, and is used to study back- 
ground from e + e~ -^qq (q — u, d, s, c) continuum events. 
The BABAR experimental program includes the measure- 
ment of the angle f3 through the measurement of time- 
dependent CP- asymmetries, A C p as discussed in Ref. [7]. 
A C p is defined as 



A C p(t) 



_ 7V(B U (t)^/)-7V(B°(<)^/) 



N(B"(t)^f)- 
= S sin(Amrft) 



N(B°(t) -> /) 
- C cos(Am d t), 



(1) 



— — 

where N(B (t) — > /) is the number of B that decay 
into the CP-eigenstate / after a time t and Amj is the 
difference between the B mass eigenstates. The sinu- 
soidal term describes interference between mixing and 
decay and the cosine term is the direct CP asymmetry. 
In Eq. 1, A(B°(t) -> J) (A{B°(t) -» /)) is the decay 
amplitude of B (B°) to the final state / (/). 

In this article, the current status of measurements of 
sin2/3 and cos2/3 from b — > ccs decays at BABAR are dis- 
cussed. All results are final unless otherwise stated. Ad- 
ditional results on CP violation measurements in B to 
charm decays at BABAR can be found in Ref. [8] . 
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The determination of j3 from b — > ccs decay modes cur- 
rently provides the most stringent constraint on the uni- 
tarity triangle. For these decay modes, the CP violation 
parameters S and C are Sb^ccs = —fjf sin2/? and C^ccs 
= 0, where 77/ is —1 for (cc)K Q s decays (e.g. J/tpK®, 
tp{2S)K^, XdK° s , r) c K Q s [9]) and r\ s is +1 for the {cc)K° L 
(e.g. J/^K° L ) state. We use the value 77/ = 0.504 ±0.033 
for the J/^jK* (K* -> K°n°) final state since it can be 
both CP even and CP odd due to the presence of even 
and odd orbital angular momentum contributions [10]. 
These modes have most recently been used to measure 
sin2/3 using a sample of 347.5 x 10 6 T(4S') ~> BB de- 
cays [11]. This result is preliminary. 

In addition to the CP modes described above, a 
large sample Bfl av of B° decays to the flavor eigenstates 
D^-h+(h+ = tt+,p+, and a+) and J/tPK*°(K*° -» 
K + it~) is used for calibrating the flavor tagging per- 
formance and At resolution. Studies are performed to 
measure apparent CP violation from unphysical sources 
using a control sample of B + mesons decaying to the 
final states J/ipK^ + , ifj(2S)K+, Xc\K + , and rj c K+ . 
The event selection and candidate reconstruction are un- 
changed from those described in Refs. [7, 12, 13] with 
the exception of a new ij c Kg event selection based on the 
Dalitz structure of the r/ c — > K®K + tt~ decay. We calcu- 
late the time interval At between the two B decays from 
the measured separation Az between the decay vertices 
of B rec and Bt ag along the collision (z) axis [7]. The z po- 
sition of the B rec vertex is determined from the charged 
daughter tracks. The B tag decay vertex is determined 
by fitting tracks not belonging to the B roc candidate to 
a common vertex, employing constraints from the beam 
spot location and the B rcc momentum [7]. Events are 
accepted if the calculated At uncertainty is less than 
2.5 ps and |At| is less than 20 ps. The fraction of events 
satisfying these requirements is 95%. 

At the T(AS) resonance, Aqp is extracted from the 
distribution of the difference of the proper decay times, 
t = tcp — ttag , where tc p refers to the decay time of the 
signal B meson (Bcp) and t tag refers to the decay time of 
the other B meson in the event (Btag)- Multivariate al- 
gorithms are used to identify signatures of B decays that 
determine ( "tag" ) the flavor of the B tag at decay to be ei- 
ther a B° or B°candidate. These algorithms account for 
correlations among different sources of flavor information 



and provide an estimate of the mistag probability for each 
event. Each event whose estimated mistag probability is 
less than 45% is assigned to one of six tagging categories. 
The Lepton category contains events with an identified 
lepton; the remaining events are divided into the Kaon I, 
Kaon II, Kaon-Pion, Pion, or Other categories based on 
the estimated mistag probability. For each category i, 
the tagging efficiency E{ and fraction Wi of events hav- 
ing the wrong tag assignment are measured from data. 
The figure of merit for tagging is the effective tagging 
efficiency Q = E 4 £ »( 1 - 2w tf = ( 30 - 4 ± °- 3 ) %, where 
the error shown is statistical only. 

With the exception of the J/ipK® mode, we 
use the beam-energy substituted mass ities = 
\A-^bcam) 2 — (Pb) 2 *° determine the composition of our 
final sample, where E£ eam and p* B are the beam energy 
and B momentum in the e + e~ center-of-mass frame. 
For the J/ipK® mode we use the difference AE be- 
tween the candidate center-of-mass energy and £J eam . 
We use events with m ES > 5.2GeV/c 2 (AE < 80MeV 
for J/ipK° L ) in order to determine the properties of the 
background contributions. We define a signal region 
5.27 < m ES < 5.29GeV/c 2 [\AE\ < lOMeV for J/ipK° L ) 
that contains CP candidate events that satisfy the tag- 
ging and vertexing requirements as listed in Table I. 

TABLE I: Number of events iVt ag in the signal region after 
tagging and vertexing requirements, signal purity P includ- 
ing the contribution from peaking background, and results of 
fitting for CP asymmetries in the Bcp sample and various 
subsamples. In addition, results on the Bna, v and charged B 
control samples test that no artificial CP asymmetry is found 
where we expect no CP violation (sin2/3 = 0). Errors are sta- 
tistical only. The signal region is 5.27 < ities < 5.29GeV/c 2 
(\AE\ < lOMeV for J/ipK°). 



Sample 




P(%) 


sin2/3 


Full CP sample 


11496 


76 


0.710 ±0.034 




6028 


92 


0.713 ±0.038 


J/ipK° L 


4323 


55 


0.716 ±0.080 


J/iPK*°{K*° -> K° s -k°) 


965 


68 


0.526 ±0.284 


1999-2002 data 


3084 


79 


0.755 ± 0.067 


2003-2004 data 


4850 


77 


0.724 ± 0.052 


2005-2006 data 


3562 


74 


0.663 ± 0.062 




4076 


96 


0.715 ± 0.044 


J^K% (JC<) -> ttV) 


988 


88 


0.581 ±0.105 


l/>{2S)K% (K° s - 7T + 7T-) 


622 


83 


0.892 ± 0.120 


XciK° s 


279 


89 


0.709 ±0.174 


r) c K% 


243 


75 


0.717 ±0.229 


Lepton category 


703 


97 


0.754 ± 0.068 


Kaon I category 


900 


93 


0.713 ± 0.066 


Kaon II category 


1437 


91 


0.711 ±0.075 


Kaon-Pion category 


1107 


89 


0.635 ±0.117 


Pion category 


1238 


91 


0.587 ±0.175 


Other category 


823 


89 


0.454 ± 0.469 


Bfl av sample 


112878 


83 


0.016 ±0.011 


B + sample 


27775 


93 


0.008 ±0.017 



2 

We determine sin2/3 with a simultaneous maximum 
likelihood fit to the At distribution of the tagged Bcp 
and -Bflav samples. There are 65 free parameters in the 
fit: sin2/3 (1), the average mistag fractions w and the 

differences Aw between B° and B mistag fractions for 
each tagging category (12), parameters for the signal At 
resolution (7), parameters for CP background time de- 
pendence (8), and the difference between B° and B re- 
construction and tagging efficiencies (7); for Bfl av back- 
ground, time dependence (3), At resolution (3), and 
mistag fractions (24). For the CP modes (except for 
J/ipK°), the apparent CP asymmetry of the non-peaking 
background in each tagging category is allowed to be a 
free parameter in the fit. We fix t b o = 1.530 ps, Am^ = 
0.507 ps" 1 [14], |A| = 1, and Ar d = 0. The fit to the B C p 
and B flav samples yields sin2/3 = 0.710±0.034±0.019 [15]. 
Figure 1 shows the At distributions and asymmetries in 

yields between events with _B°tags and B tags for the 
r]f = —1 and t]f = ±1 samples as a function of At 
overlaid with the projection of the likelihood fit result. 



~400 




FIG. 1: a) Number of = -1 candidates (J/ipK%, ip(2S)K° s , 
XciK®, and r/ c Kg) in the signal region with a £>°tag (N B o) 

and with a B tag (N-^o), and b) the raw asymmetry (N B o — 
N-go)/(N B o + N-go), as functions of At. Figures c) and d) 
are the corresponding distributions for the Tjf = ±1 mode 
J/tpK°. All distributions exclude Other-tagged events. The 
solid (dashed) curves represent the fit projections in Atfor 
B°(B ) tags. The shaded regions represent the estimated 
background contributions. 

We perform a separate fit with only the cleanest f]f = 
— 1 sample, in which we treat both |A| and sin2/3 as 
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free parameters. We do not use the modes J/tpK*° 
and J/ipK® to minimize the dependence of the results 
on the background parametrization. We obtain |A| = 
0.932 ± 0.026 ± 0.017. The updated value of sin2/3 is 
consistent with the current world average [16] and the 
theoretical estimates of the magnitudes of CKM matrix 
elements in the context of the SM [17]. 

Ref. [18] presents a model-independent study of this 
shift using the measurements of B° —> J/iJ;ir [19] to 
quantify the effect of contributions from penguin oper- 
ators and long-distance contributions from penguin con- 
tractions. They find that the deviation of the measured 



Sep term from sin2/3, A5 



S 



sin2/3 



0.000 ± 0.017 which is comparable to the systematic er- 
ror from our previous publication [12]. The theoretical 
estimates ofASj^ K o arc 0(1Q- 3 ) [20] and 0(1Q- 4 ) [21]. 



III. cos2/3 FROM b -> ccs DECAYS. 

The analysis of b — > ccs decay modes imposes a con- 
straint on sin2/? only, leading to a four-fold ambiguity in 
the determination of f3. This ambiguity can leave pos- 
sible new physics undetected even with very high preci- 
sion measurements of sin2/3. Additional constraints are 
obtained from the ambiguity-frcc measurement of cos2/3 
using the angular and time-dependent asymmetry in 
B° — > J/ip K* decays and the time-dependent Dalitz plot 
analyses of B° ->' D*°h° and B° D*+D*~K°. The 
BABAR B° D*°h° analysis is described in Refs. [22] 
and [23]. 



D* + D*-K% 



Both B° and B" mesons can decay to the same final 
state D* + D*~K° via a process dominated by the sin- 
gle weak phase VT-emission b — ► ccs transition. Possible 
penguin contributions are neglected and it is therefore 
assumed that there is no direct CP violation. According 
to Rcf. [24] , the decay can be divided into two half Dalitz 
planes s + < s~ and s + > s~ , where s + = m 2 (D* + K®) 
and s~ = m 2 (D*~ K®), such that the time-dependent 
decay rate asymmetry of B° 



D*+D*-K° s is 



A(t) 



F , r = Vy-r cos(Am d t) - 

i-gO+L B a ,J a 

' sin2/3 + r] y ^ s2 cos2/3 ) sin(Ara<ii), 



Jo 



Jo 



where i] y = — 1(+1) for s + < s~(s + > s~). The pa- 
rameters Jo,J c ,J s i and J s2 are the integrals over the 
half Dalitz phase space with s + < s~ of the functions 
|a| 2 + |a| 2 , |a| 2 — |a| 2 , Tle(aa*) and Xm(aa*), where a 
and a are the decay amplitudes of B° -» D*+ D*~ K° s 

and B° — > D* + D*~ K®, respectively. If the decay 
B° — > D* + D*~ K® has only a non-resonant component, 



the parameters J S 2 — and J c are at the few percent 
level [24]. The CP asymmetry can be extracted by fit- 
ting the B° time-dependent decay distribution. The mea- 
sured CP asymmetry is sin2/3 multiplied by a factor of 
2J s i/J because the final state is an admixture of CP 
eigenstates with different CP parities. In this case, the 
value of the dilution factor 2J s i/J is estimated to be 
large [24], similar to the decay B Q -> D*+D*~. The 
situation is more complicated if intermediate resonances 
such as Dfj are present. In this case, the parameter 
J S 2 is non-zero and J c can be large. The resonant com- 
ponents are expected to be dominated by two P-wave 
excited D sl states [24]. One such state is £^(2536) that 
has a narrow width and does not contribute much to J S 2 ■ 
The other resonant state is predicted in the quark 
model [25] to have a mass above the D* + K® mass thresh- 
old with a large width. In this case, the J S 2 can be large. 
Therefore by studying the time-dependent asymmetry of 
B° -> D*+D*-K° in two different Dalitz regions, the 
sign of cos 2/3 can be determined for a sufficiently large 
data set using the method described in Refs. [24, 26]. 
This would allow the resolution of the (i — > n/2 — [3 
ambiguity despite the large theoretical uncertainty of 
Usij Jo- However, one of the expected P-wave may 
be the newly discovered D+ (2317) or £>+ 7 (2460). These 
states are below the D* + K® mass threshold, so they will 
not contribute to the decay BP_ -> D*+D*-R°. The 
fits to 208.7 x 10 6 r(45) -» BB decays yield: J c /J = 
0.76 ± 0.18 ± 0.07, (2J sl /J ) sin2/3 = 0.10 ± 0.24 ± 0.06, 
and (2J s2 /J ) cos2/3 = 0.38 ± 0.24 ± 0.05 [27]. The mea- 
sured value of J c / Jo is significantly different from zero, 
which, according to Ref. [24], may indicate that there 
is a sizable broad resonant contribution to the decay 
B° —> D* + D*~ K® from an unknown state with an 
unknown width. Under this assumption then the mea- 
sured value of 2J S 2/Jo implies that the sign of cos2/3 is 
preferred to be positive at a 94% confidence level. 



B. B° -> J/^K ,0 (K*° -> K° s tt°). 

BABAR has measured the sign of cos2/3 in a time- 
dependent angular analysis of 104 B a — > J/ipK*°(K* — > 
K °tt°) decays in 88.0 x 10 6 Y(4S) -> BB of data recorded 
between 1999 and 2002 [28]. Interference between decays 
to CP-even (L=0,2) and CP-odd (L=l) final states give 
terms proportional to cos2/3 in the decay rate. Strong 
phase differences and transversity amplitudes A, that ap- 
pear also in these terms, have been separately measured 
in a time-integrated angular analysis of B^ — ► J/ipK*^ 1 
and J/tP K*°(K* ^ K+tt-) decays: 



6\\ - S = 

5±-6 = 

\A\ 2 = 

l^lll 2 = 

\A i I 2 = 



(-2.73 ± 0.10 ± 0.05) rad, 
(+2.96 ± 0.07 ± 0.05) rad, 
0.566 ±0.012 ± 0.005, 
0.204 ±0.015 ± 0.005, 
0.230 ±0.015 ±0.004. 



4 



The analysis in principle allows a second solution for the 
strong phase differences, leading to a sign ambiguity in 
cos2/3. This ambiguity has been resolved with the inclu- 
sion of S-wave Kir final states in the analysis. The in- 
terference between the S-wave and P-wave contributions 
gives additional terms in the decay rates with a clear de- 
pendence on the Kir mass due to the resonance shapes. 
The other solution for the strong phase differences can 
be excluded as leading to an unphysical dependence of 
the strong phase differences on the Kir mass [29]. Us- 
ing the values from Eq. 2, and fixing sin2/3 to 0.731 1 , 
the fit to the B° -> J/ijK* a (K*° -> K° s 7r°) sample gives 
cos2/3 = +2.72^°;™ ± 0.27. By comparing this result 
with the outcome of fits to 2000 data-sized Monte Carlo 
samples, the sign of cos 2/3 is determined to be positive 
at the 86% confidence level, in agreement with Standard 
Model expectations. 



IV. CONCLUSIONS. 

When the BABAR measurement of sin2/3 using b —> ccs 
decays is combined with the most recent Belle result 

1 This was the world average at the time - see Rcf. [30] . 



described in Ref. [31] then the world average value of 
sin2/3 from b -> ccs decays is sin2/3 = 0.674 ± 0.026. 
The combined constraint on (3 in the ~p-rj plane from 
the BaBAR and Belle b —> ccs charmonium + K° me- 
son analyses [11, 31], the B a -> J/^jK* [28, 32], B° -> 
D*°h° [22, 33] and B° -> D*+D*-K° s [27] analyses 
strongly favour the solution (3 =21.1±1.0° where cos2/3 
is positive [16]. 
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